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The spontaneous formation of droplets via dewetting of a thin fluid film from a solid substrate 
allows for materials nanostructuring, under appropriate experimental control. While thermal fluc¬ 
tuations are expected to play a role in this process, their relevance has remained poorly understood, 
particularly during the nonlinear stages of evolution. Within a stochastic lubrication framework, we 
show that thermal noise speeds up and substantially influences the formation and evolution of the 
droplet arrangement. As compared with their deterministic counterparts, for a fixed spatial domain, 
stochastic systems feature a smaller number of droplets, with a larger variability in sizes and space 
distribution. Finally, we discuss the influence of stochasticity on droplet coarsening for very long 
times. 

PACS numbers: 47.61.-k, 05.10.Gg, 68.08.Bc, 68.15.+e 


The interplay between stochastic fluctuations and non¬ 
linear interactions can induce highly nontrivial effects in 
spatially extended systems [1] at the nanoscale [2]. For 
instance, noise can rectify the direction of material trans¬ 
port, as for diffusing particles under asymmetric forces 
[3]. When a characteristic pattern emerges from a ho¬ 
mogenous state [4] , fluctuations can even enhance (rather 
than hinder) spatial order, or modify the rate at which 
typical pattern sizes increase with time (coarsening), as 
e.g. for evolving atomic steps at epitaxial surfaces [HE]. 


A natural context in which fluctuations are expected to 
be relevant is nanoscale fluid flow. Although the contin¬ 
uum framework provided by the Navier-Stokes equations 
is physically valid down to surprisingly small scales (c^ 1 
nm) w, the atomistic nature of the fluid medium is ex¬ 
pected to play an increasingly important role as physical 
scales are reduced. Note that in this process the surface- 
to-volume ratio also becomes ever more favorable [9], so 
that free surface flows m provide conspicuous instances 
for noise effects. Thus, a number of interfacial processes 
have been seen to depend critically on the occurrence of 
fluctuations, such as e.g. the breaking of nanojets m 
ng or of liquid threads m- In addition, it is known that 
thermal noise in the fluid velocity field changes the value 
of the contact angle under partial wetting conditions [15] , 
and enhances the spreading of droplets in surface-tension 
m and gravity m driven systems, as well as the rup¬ 
ture of thin dewetting films iniiiH]. 


Indeed, in dewetting experiments carried out using 
polymer m or liquid metal films m, dynamics and mor¬ 
phologies are not reproducible by deterministic frame¬ 
works. Early rupture times and irregular patterns sug¬ 
gest that fluctuations play a strong dynamical role. 
Working in the long wave (lubrication) approximation 
to free surface flow [ 21 ], previous works improving de¬ 
terministic models [TBHlQ] [ 22 ] cast the time evolution of 
the system in the form of a stochastic evolution equation 


(SEE) for the thickness h of the thin fluid film, namely. 


rjdth = dx 



n(/i)] + a e(x, t) I , 
( 1 ) 


where r] is viscosity, 7 is surface tension, e(x,t) is a 
Gaussian white noi se of zero mean and unit variance, 
and (j = y^r]kBT/3 depends on temperature T. In Eq. 

n(/i) = —d^(h)/dh is the disjoining pressure that 
accounts for fluid-solid interaction, with 4>(/i) the in¬ 
terface potential [23]. A power law is commonly used, 
n(/i) = K[(h^/h)^ — (h^/h)^], where n is proportional to 
the Hamaker constant and is the precursor film thick¬ 
ness isaEii, corresponding to the minimum of the po¬ 
tential. The equilibrium contact angle emerges from the 
interplay between surface tension and disjoining pressure, 
see e.g. [25] . 

As shown in CHI EH], for short-times Eq. 0 predicts a 
morphological instability [4]: surface tension, in compe¬ 
tition with the destabilizing disjoining pressure, selects a 
typical length scale. A, of surface undulations. In a pro¬ 
cess reminiscent of domain coarsening in phase separation 
or spinodal decomposition of binary mixtures [28], and 
as seen in the experiments m, this scale increases non- 
trivially with time in a form which cannot be accounted 
for deterministically [i.e., for cr = 0 in Eq. Q]. However, 
the study in m was limited to the linear regime, where 
perturbations of the flat film are small compared to its 
thickness and the morphology is dominated by capillary¬ 
like surface modes. Pressing open questions then are 
if and how do the nonlinearities, responsible for actual 
droplet formation in Eq. 0 > modify this time evolution. 
Moreover, at very long times deterministic droplets are 
known to undergo coarsening into a single-drop morphol¬ 
ogy ESI EH- Hence, one can ponder whether thermal 
fluctuations modify the coarsening behavior in this case. 
Eor instance, noise is known to modify the coarsening law 
of the ID Cahn-Hilliard equation, a paradigmatic model 
in the context of spinodal decomposition EH] [29] . 

In this Letter we study the effect of thermal fluctua- 
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tions on the formation and evolution of droplets under 
partial wetting conditions. To this end, we study nu¬ 
merically the SLE [Eq. 0 ] in the nonlinear regime. We 
find that thermal fluctuations unambiguously speed up 
the nonlinear process of droplet formation. Moreover, 
we show that, as compared to the deterministic case, 
noise increases heterogeneity in droplet sizes and posi¬ 
tions, while seemingly not affecting the coarsening pro¬ 
cess expected for asymptotically long times PHUT]. 

To study the behavior of the solutions of Eq. 0 in 
nonlinear regime, we have carried out large-scale numer¬ 
ical simulations using a scheme [30] that ensures non- 
negativeness of the solution for all times, if the initial 
film is positive in the full domain m- Specifically, 
our algorithm is based on the standard implicit (Crank- 
Nicholson) discretization [30|, where surface tension is 
treated implicitly, while n(h) is treated explicitly; we 
employ zero-flux boundary conditions. The stochastic 
term in Eq. 0 is also dealt with explicitly, within the 
Stratonovich interpretation m- 

In our simulations we consider a nondimensional ver¬ 
sion of Eq. Q obtained by defining h = hjhc, x = xjhc^ 
and t = t/tc, where he is a typical film thickness and 
tc = 3r]hc/^ [32]. This leads to non-dimensional ampli¬ 
tudes a = and to ^ = Khe/^\ we use the 

exponents (n^m) = (3,2) as in e.g. [25]. We perform the 
deterministic and stochastic simulations of Eq. 0 using 
a precursor thickness = 0.01 and the same random 
initial condition^ namely, random values of the thick¬ 
ness with non-dimensional average Hq = 0.1 and vari¬ 
ance 10~^hQ. The contact angle is set to 50° in the ex¬ 
pression k = 2(1 — cos^)//i* [25], leading to = 72; 
within the long wave theory implementation, the actual 
contact angle (measured by the slope of the tangent line 
passed through the drop profile through the point of in¬ 
flection) is smaller and is close to 25°. The spatial step 
size dx = h^; this choice is known to lead to accurate 
results [30]. The temporal step size is adaptive, follow¬ 
ing the approach described in m- This allows to obtain 
converged results with reasonable computational effort; 
the use of such adaptive time stepping is particularly im¬ 
portant for the purpose of carrying out simulations for 
long times where coarsening effects become relevant. We 
use a large domain size L ^ 31 A, where A is the most 
unstable wavelength obtained by linear stability analy¬ 
sis of deterministic version of Eq. Q, discussed further 
below. Combined with a large number of realizations, 

200 , such a domain allows to obtain statistically mean¬ 
ingful results. The specified parameter values are closely 
related to the polymer films studied in [19], where the 
characteristic film thickness is 4 nm (he = 40 nm, so that 
ho = 0.1), while 7 = 0.03 N/m and the Hamaker con¬ 
stant ^4 = 2 • 10“^^ J yield a contact angle in the 15 — 20° 
range. Eor these parameter values the non-dimensional 
noise strength a 10“^ corresponds to T = 50 — 60° C. 
On the other hand, for the liquid metal thin films consid¬ 
ered in [ 20 ], 7 = 1.3 N/m, T = 2000 K, and he G [50,150] 




FIG. 1. Space-time plot of droplet formation and evolution 
as predicted by Eq. for a = 0 (a) and a — 10“^ (b), for 
the same parameter values and initial conditions, see main 
text. Brighter (darker) color corresponds to larger (smaller) 
values of the film thickness h{x,t). 


nm, leading to d G [10“^/^, 10“^]. 

Erom now on, and unless otherwise stated, we work in 
dimensionless units and remove hats for notational sim¬ 
plicity. Eigure shows examples of the time evolution 
predicted by Eq. Q in the deterministic and stochas¬ 
tic cases. Well-defined droplets (clear bands) emerge af¬ 
ter a rupture time of roughly tr,det = 180 (^r,sto = 80) 
time units in the deterministic (stochastic) system. In 
the latter case there is a substantial amount of droplet 
merging around that time, after which activity decreases. 
Comparing both panels, we immediately observe that the 
width of the droplets (clear bands) and their locations 
are much more regular in the deterministic than in the 
stochastic case. 

Although some spatial modulation can be seen for ear¬ 
lier times in Eig. the system behavior is less visually 
clear. However, at such times one can resort to linear 
stability analysis [HlTlllS]. The time evolution of the 
system is conveniently described by the structure factor 
Sq = (|hg(t)p), which within linear approximation can 
be analytically obtained [13 Ha], 


5 , 


(27r)^ 


So{q)e 


2(jo{q)t 


+ 


2 uj{q) \ 



( 2 ) 

Here, hq{t) is the Eourier cosine transform [33] of h{x,t) 
for wavenumber g, So{q) is the initial structure factor, 
ho is a film thickness, and the growth rate is given by 
the dispersion relation uj{q) = h^q^ (2gQ — g^) /3. Here 
go = —n'(ho)/2. The wavelengths of unstable perturba¬ 
tions correspond to g G [0,\/2go], for which uj{q) > 0. 
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FIG. 2. Time evolution of the position of the main maximum, 
Qm, of the structure factor. The dashed black line indicates 
the deterministic linear prediction, qo. Black circles (squares) 
correspond to predictions from Eq. 0 for cr = 10“^ (a = 0). 
Red up [blue down] triangles provide the position of gm,sto(t) 
[^m,det(^)] as obtained in numerical simulations of Eq. M for 
a = 10“^ (cr = 0). Rupture times are signalled by arrows. All 
results are obtained by averaging over 200 noise realizations. 
Inset: Number of droplets for different noise amplitudes at 
t = 220. All lines are guides to the eye. 


Starting from an initial condition with mean ho^ the de¬ 
terministic system very quickly selects the wavenumber 
^m,det = ^0 for which the growth rate uj{q) reaches its 
positive maximum, see black squares and blue triangles 
in Fig. where we plot the time evolution of the value 
of wavenumber q^ at which the main maximum of Sq 
occurs. For our parameter choice, qo = 2.464. Within 
linear approximation, this sets the length scale of the 
pattern, A = 27r/go = 2.546, namely, the average size of 
surface undulations. In contrast, stochastic systems ini¬ 
tially develop nontrivial short lengthscale (large q) corre¬ 
lations, so that Sq{t) displays a maximum for a wavenum¬ 
ber value gm,sto which decreases with time towards the 
deterministic value qo, see Fig. and [32]. This is the 
process described in m as coarsening. Note that, as 
mentioned above, droplets have not yet formed; as seen in 
[32| , for these times the film morphology remains largely 
a small-amplitude sinusoid. In addition, for stochastic 
simulations, ^m,sto > Qo] as we will see, this inequality 
does not hold in the nonlinear regime. 

Within the range of validity of the linear approxima¬ 
tion, the film develops unstable modes that remain in¬ 
dependent of one another. If the linear predictions were 
applicable to long times, then the number of drops even¬ 
tually formed would be essentially fixed by the linear 
value A = 2.546, since Sq is characterized by a well 
defined peak around q = Qm^ see black lines in Fig. 
§ However, experiments show that the distribu¬ 
tion of droplet sizes is relatively wide. Droplet forms 
differ strongly from smooth sinusoids, and they interact 
non-trivially (e.g. through merging and coalescing) dur- 




FIG. 3. (a), (c): Surface morphologies from simulations of 

Eq. Q for cr = 0 (blue line) and a = 10“^ (single realization, 
thick red line) at t = 120 (a) and 220 (c). (b), (d): Structure 
factor averaged over 200 noise realizations, at t = 120 (b) 
and 220 (d) for a = 0 (blue squares) and a = 10“^ (red 
circles). The thick black lines in (b) provide the corresponding 
analytical predictions from Eq. Q. Thin lines are guides to 
the eye. 


ing their evolution. On long time scales, the number of 
drops needs to reduce, as a single larger drop is ener¬ 
getically more favorable than two smaller ones [261127]. 
Actually, the most stable configuration of the system is 
a single droplet, since the evolution described by Eq. 0 
drives the system to the minimum of the interface Hamil¬ 
tonian l-Llh] = f dx [4>(/i) + 7(9a,/i)^/2] [TS] , 

Hence, we next need to address droplet formation for 
times t > 60, see Fig. 2, away from the linear regime. As 
seen in the animation provided at [32], nonlinear effects 
indeed set in for t 60. Thus, the deterministic structure 
factor develops higher harmonics, while the stochastic Sq 
also departs clearly from the linear solution, Eq. Q, see 
Fig-i for two sample times. The higher harmonics are 
at least one order of magnitude smaller than the main 
peak [32] , so that they barely influence later stages of the 
evolution. In addition, the rupture time at which well- 
defined droplets form is much shorter for the stochastic 
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(^r,sto — 80) than in the deterministic (tr^det — 180) case, 
see [32] and also Fig.j^a) for t = 120, where droplets have 
appeared in the former case, but not yet in the latter. 

After rupture, the Sq distribution broadens around the 
main peak both in the stochastic and in the determinis¬ 
tic systems, and for values of q on both sides of see 
Fig. [lb ,d) and [32|. Moreover, there is an additional 
boost in the rupture process so that stochastic droplets 
create faster than one would expect using the linear the¬ 
ory: Note that rupture times are signaled by a kink in the 
corresponding qm{t) data. At rupture, nonlinear ripen¬ 
ing of droplets takes place, namely, a decrease of q^ with 
time, which is more pronounced and occurs earlier in the 
stochastic system. In contrast to linear predictions, the 
deterministic system also undergoes a similar, albeit de¬ 
layed process. We conjecture that disorder in droplet 
positions favors merging of nearby drops, inducing more 
rapid decrease of q^ in the stochastic system. 

Also, for any cr > 0 , once the drops are well formed the 
decrease in qni{t) with time slows down. On average, the 
value of q^ which is eventually achieved (say, for t > 220 ) 
is smaller for a 7 ^ 0. This behavior implies a smaller num¬ 
ber of drops for a fixed domain in the stochastic system, 
see e.g. Fig. [^c). Also recall Fig. where substantial 
drop merging is seen for this case during times from rup¬ 
ture up to t 120. We note that, for the time scales 
considered so far, the final number of droplets decreases 
when the noise intensity (say, temperature) increases, see 
the inset of Fig. 

Figure showing the distribution of drop heights and 
their distances, illustrates a further significant difference 
between stochastic and deterministic evolution: stochas- 
ticity indeed leads to much wider droplet distributions, 
and therefore to much more irregular patterns. Two sam¬ 
ple morphologies are compared in Fig.j^c). Also, the in¬ 
set in Fig.j^b) shows that the width at half maximum of 
the Sq distribution. A, is an increasing function of noise 
amplitude, as expected. This finding may be of signifi¬ 
cant importance in applications, where regularity of the 
distribution of drops is often desired. Our results sug¬ 
gest that decreasing noise amplitude may be the key to 
achieve this goal. 

Finally, we consider much larger time scales, to probe 
the convergence of the stochastic system to the equilib¬ 
rium single-drop solution. Up to the times discussed so 
far, the decrease of qm{t) seems mostly induced by droplet 
coalescence. This introduces relatively large distances 
among remaining units, recall Fig. for long times. 
For still longer times, droplet interaction occurs mostly 
through the precursor layer, inducing non-interrupted 
coarsening of the pattern into a single drop morphology 
[26|. For well-separated droplets and a = 0, analytical 
predictions actually exist for the decrease of the number 
of droplets N{t) with time [^. We have considered the 
evolution predicted by Eq. Q for a 7 ^ 0 at very long 
times, up to t = 10 ^. Computational feasibility requires 
a larger precursor thickness h* = 0.04 and smaller noise, 
a = 10“^/^. Our results indicate that fluctuations do 
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FIG. 4. Distribution functions of drop heights (a) and inter¬ 
drop distances (b) at time t = 220 for cr = 0 (blue squares) 
and cr = 10~^ (red circles). Inset: width of the main peak of 
Sq dX t — 220, as a function of noise amplitude. All lines are 
guides to the eye. 



shorten significantly the time scales on which coarsening 
occurs. However, they become less relevant with increas¬ 
ing time, to the extent that the asymptotic behavior of 
N{t) is not modified with respect to the deterministic 
case, at least within the accuracy of the results, see [34] . 
Thus, droplet coarsening counts among phenomena for 
which noise does not change the coarsening universality 
class [28] of the corresponding deterministic system. 

In summary, we have shown that stochastic effects due 
to thermal noise may play a significant role in dewetting 
of thin fluid films, in each of the three stages of evolu¬ 
tion considered. For very early times, stochasticity leads 
to a decrease of the most unstable wavenumber, g'm,sto 5 
down from the values that are large compared to the 
deterministic one, go; however, within this stage gm,sto 
remains larger than go- After this, noise triggers an ear¬ 
lier onset of nonlinear effects, inducing a shorter rupture 
time. At these time scales, stochasticity leads to droplet 
coarsening, in the sense that gm,sto < contrast to 

the linear regime. Finally, for much longer times, fluc¬ 
tuations speed up the coarsening process that will ulti¬ 
mately lead to formation of the single-drop, energetically 
favored state. Qualitatively, the deterministic coarsening 
law for the number of drops remains unchanged. How¬ 
ever, quantitatively the time scales involved in this long¬ 
time coarsening process are significantly influenced by 
noise, and we conjecture that stochastic effects may be 
observable in careful experiments carried out with fluid 
films of nanoscale thickness. 
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Supplemental Material 

In this Supplemental Material we provide results from numerical simulations of the stochastic thin film equation 
discussed in the main text, which in the corresponding (hatted) dimensionless units reads 

dth = V-{h^W[-V'^h-n{h)]+ah^/'^€{x,t)^, (3) 

where e is an uncorrelated Gaussian noise with zero mean and unit variance, and the hats have been removed to 
simplify the notation. 

The film morphology h(x,t), the structure factor 5'g(t), together with the histograms of individual drop heights 
and inter-drop distances, are all shown in Movie SI for times t G [1,250], which includes deterministic (a = 0) and 
stochastic (a = 10“^) simulations of Eq. These simulations have been performed using as initial condition a 
film of non-dimensional height ho = 0.1 (perturbed by Gaussian white noise with 10“^ amplitude), with a precursor 
thickness =0.01. Other parameters are described in the main text. 

In order to probe the long-time evolution of the droplet pattern, in which coarsening takes place mainly through the 
precursor film that communicates individual droplets m, we have considered a larger precursor thickness, = 0.04. 
Gomputational feasibility has also required us to set a = 10“^/^. The time evolution of the film morphology under 
these conditions can be assessed in Movie S2, which corresponds to stochastic simulations. Figure is a plot of the 
number of droplets vs time, N{t), in this long-time coarsening regime, averaged over 40 realizations of the noise. 
While the decay of N{t) may not be far from analytical predictions and simulations for very large discrete models 
[26] . our simulations seem to support the irrelevance of noise at sufficiently long times t > 2000. 



FIG. 5. Log-log plot showing the time evolution of the number of droplets vs time, for long times up to t = 10000. Blue squares 
(red circles) correspond to deterministic {a = 0) and stochastic {a = simulations, respectively. For reference, the dashed 

lines correspond to the power-law decay N{t) ~ found for very large deterministic systems in [26]. Solid blue and red 

lines are guides to the eye. 









